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ABSTRACT 

The aperture mass statistic is a common tool used in weak lensing studies. By convolv- 
ing lensing maps with a filter function of a specific scale, chosen to be larger than the 
scale on which the noise is dominant, the lensing signal may be boosted with respect 
to the noise. This allows for detection of structures at increased fidelity. Furthermore, 
higher-order statistics of the aperture mass (such as its skewness or kurtosis) , or count- 
ing of the peaks seen in the resulting aperture mass maps, provide a convenient and 
effective method to constrain the cosmological parameters. In this paper, we more 
fully explore the formalism underlying the aperture mass statistic. We demonstrate 
that the aperture mass statistic is formally identical to a wavelet transform at a spe- 
cific scale. Further, we show that the filter functions most frequently used in aperture 
mass studies are not ideal, being non-local in both real and Fourier space. In contrast, 
the wavelet formalism offers a number of wavelet functions that are localized both 
in real and Fourier space, yet similar to the optimal aperture mass filters commonly 
adopted. Additionally, for a number of wavelet functions, such as the starlet wavelet, 
very fast algorithms exist to compute the wavelet transform. This offers significant 
advantages over the usual aperture mass algorithm when it comes to image processing 
time, demonstrating speed-up factors of ~ 5 — 1200 for aperture radii in the range 2 
to 64 pixels on an image of 1024 x 1024 pixels. 

Key words: gravitational lensing: weak - methods: data analysis - cosmology: cos- 
mological parameters - cosmology: dark matter 



1 INTRODUCTION 



tool in weak gravitationa l lensing s t udies , being useful for 
detec tion of structures (I Schneider! 1 19961 : ISchneider et al.l 
1998) and for the determination of cosmological parameters 
through related statistics such as its varian ce, skewness and 



kurto si s, or throug h peak s tatistics (e.g. Schneider et al.l 
19981: IJarvis et all 120041 : iKilbinger fc Schneider! 120051 : 
Dietrich fc Hartlad 120101 ). Formally, the aperture mass 
technique consists of applying a filter, defined on circular 
apertures, to a map of measured shears. 

The strength of the aperture mass technique lies in the 
fact that the structures responsible for the lensing signal 
are dominant over the noise at certain scales. Choosing an 
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aperture filter whose shape traces that of the expected signal 
and with a scale similar to that of the expected structures 
results in an aperture mass map with optimal signal to noise. 
Moreover, assessment of the noise properties of a given M ap 
reconstruction is straightforward, carried out by randomis- 
ing the input data a large number of times and computing 
the variance of the resulting M ap reconstructions, <tm ■ The 
signal-to-noise of the aperture mass map can thereby be di- 
rectly computed from the data : S = ^ p . 

However, the signal-to-noise will only be maximum in 
the aperture mass map if the filter function is tuned to fol- 
low the shape of the expected signal, which is not necessar- 
ily known a priori, and when the aperture size matches the 
scale of the structures we aim to detect. Choosing an ideal 
filter function is, consequently, somewhat challenging. Per- 
haps the largest difficulty with this method, though, is that 
evaluating the M ap statistic and its associated noise map is 



© 2009 RAS 



2 A. Leonard, S. Pires, J.-L. Starch 



time consuming, particularly for larger apertures, and this 
limits the ability to consider reconstructions of large fields, 
or across a range of scales. 

An alternative method for analysing weak lensing data 
at multiple differ ent scales simulta neously is to use the 
wavelet formalism (Star ck et al hood ). This method involves 
first converting the shear map, evaluated via measurement of 
galaxy shapes, into a convergence map in Fourier space, then 
applying a wavelet transform. One advantage of this method 
is that there exist fast algorithms to compute the wavelet 
transform. These algorithms compute simultaneously sev- 
eral wavelet bands analogous to aperture mass maps at 
dyadic scales. 

In this paper, we demonstrate that the aperture mass 
statistic is formally equivalent to the wavelet transform 
of a map of the weak lensing convergence k evaluated 
at a particular scale, for appropriate choice of filter func- 
tion. Furthermore, we demonstrate t hat the most frequently 
used aperture mass filter functions ( Schneider et al.l 1 19981 ; 
Ivan Waerbekel Il998l ; iJarvis et all l2004h do not meet all of 
the desired requirements of aperture mass filters, and are 
therefore not ideal. We propose an alternative filter function 
- the starlet wavelet - and demonstrate that this filter func- 
tion meets all the desired requirements for aperture mass 
filter functions. Moreover, we show that this function is very 
similar in form to the aperture mass filter func tion found to 
be opt imal for weak lensing skewness studies in lZhang et al.l 
(2003), and is therefore a perfect choice for lensing studies. 
Lastly, we demonstrate that the a trous wavelet transform 
algorithm , publicly available as part of the MRLens software 
packagcQ ( Starck et alhOOd ). is substantially faster than the 
aperture mass algorithm, offering up to three orders of mag- 
nitude speed-up in processing time for an 1024 x 1024 pixel 
image. 

This paper is organised as follows. In §[51 we describe the 
formalism underlying the aperture mass method, and out- 
line the desired characteristics of the associated filter func- 
tions. In § |31 we define the wavelet transform of an image, 
and demonstrate that the wavelet transform of a map of 
the weak lensing convergence k is identical to the aperture 
mass statistic evaluated at corresponding scales. In § 3] we 
describe the two most commonly-used aperture mass filter 
functions, and compare these to the starlet wavelet func- 
tion, which we show to be more ideally-suited for aperture 
mass studies. In § [5] we compare the processing time for 
the aperture mass and wavelet transform algorithms, and 
we conclude with a summary of our results in § [5] 



M a p definition 



2 THE APERTURE MASS STATISTIC 
2.1 

The aperture mass statistic is defined l|Schneiderl 1 199d ) as 
the convolution of the convergence, k, with a radially sym- 
metrical filter function, 



M ap (0) = / cftf k(0)C/(|0|) 



(1) 



By considering the relationship between the shear, 7, and 
the convergence, one can reformulate equation {TJ in terms 
of the measured shear as 



Af ap (0) = / <?0 TtWQQW 



(2) 



where 7t(i9) is the tangential component of the shear at 
position i9 relative to the centre of the aperture, Q(|i9|) is a 
second radially-symmetric function, related to f7(|i?|) by: 
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and U (0) is subject to the condition that 

r-6 

■& U(f) dl& = , 



(3) 



(4) 



with 8 being the radius of the aperture. Furthermore, Q(_*) 
and U (0) are required to go to zero smoothly at ■& — 9. It is 
also preferable that the power spectrum of f7(#) is local in 
the frequency domain, and shows no oscillatory behaviour. 
This ensures that the filter function acts as a band-pass 
filter, allowing detection of structures at the scale of interest 
only. 

The filter functions Q($) and [/($) may take any form, 
subject to the conditions outlined above, namely: 

(i) The filter function should be compensated within the 
aperture being considered. 

(ii) The filter function should be local in real space; i.e. 
it should go to zero smoothly at a finite radius, and be zero 
outside this radius. 

(iii) The filter function should be local in Fourier space, 
with no oscillatory behaviour in the power spectrum. 



ISchneiderl (| 19961 ) showed that optimal signal to noise 
will be achieved if the function Q(i9) matches the expected 
profile of jt as far as possible. This is generally not known a 
priori, so a number of generic functio ns, or families of func- 
tions, have been proposed (see, e.g 



Ivan Waerbeklll998l : [jarvis et alj|2004 ). These popular filter 



Schneider etaH Il998l ; 



http: / /irfu. cea.fr/Ast/mrlens_softwarc.php 



functions will be discussed in detail in § 0] below, where we 
will demonstrate that none of them are able to simultane- 
ously meet all three of the above criteria. 

We now consider the wavelet transform and its applica- 
tion in weak lensing. Specifically, we will demonstrate below 
that the wavelet transform at a given scale is formally equiv- 
alent to the aperture mass statistic described above. 
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2.2 Associated Statistics of M ap 

The aperture mass statistic is useful not only for generating 
filtered maps of the lensing convergence, and thereby de- 
tecting structures as a function of aperture scale, but also 
for constraining cosmological parameters through associated 
statistics of the aperture mass M ap . 

Most commonly, the variance of the aperture mass 
statistic, (A'/ap) (#)> is considered as a function of aperture 
radius 6. This measure is relat ed to the power spectrum 
P K (k) via (|Schneider et al.ll 19981 ): 

(A4) (9) = i- J d 2 kP(k)W(k9) , (5) 



where 



W{k6) = [U{k)f 



(6) 



is the power spectrum of the filter t/(i?). 

Equation ((5J can equivalently be expressed in terms of 
the two point correlat i on function of the sh e ar (see, e.g. 



Crittenden et alj |2002| ; ISchneider et all |2002| ; Ijarvis et ail 
2004 ). and this statistic helps to constrain the amplitude 



of the matter power spectrum as- 

In order to probe non-gaussianities in our cosmolog- 
ical model, higher -order statistics need to be considered. 
Ijarvis et al.1 |2004 ) present a derivation of the skewness of 
the aperture mass statistic (M^ p ) in terms of the Bispec- 
trum or, equivalently, the three-point correlation function 
of the shear. This statistic, as well as the kurtosis of the 
convergenc e and apert u re ma s s, and peak counting are con- 
sidered in iPires et aL I (|2009l . |2012l ). where the efficacy of 
these probes to discriminate between competing cosmologi- 
cal models along the Qm — as degeneracy are considered. 



3 WAVELET TRANSFORM 

The wavelet transform is a multiscale transform, where the 
wavelet coefficients of an image are computed at each posi- 
tion in the image at various different scales simultaneously. 
In one dimension, the wavelet coefficient of a f unction f(x), 
evaluated at position b and sca le a is defined as (|Starck et al.l 
1 19981 ; IStarck fc Murtagbll2006l ) : 



W(a,b) 



dx 



where ip(x) is the analysing wavelet. The analysis is analo- 
gous in 2 dimensions, with tp(x,y) — tp(x)ip(y). 

By definition, wavelets are compensated functions; 
i.e. the wavelet function tp(x) is constrained such 
that J Bl ip{x)dx = and hence, by extension 



4>(x, y)dx dy — 0. 



(8) 



According to the definition in equation 0, the continuous 
wavelet transform of an image is therefore nothing more than 
the convolution of that image with compensated filter func- 
tions of various characteristic scales. If the image f(x, y) is 



taken to be the convergence n(x,y), then for an appropri- 
ate choice of (radially-symmetric, local) wavelet, the wavelet 
transform is formally identical to the aperture mass statistic 
at the corresponding scales. 

In practice in application, we use the starlet trans- 
form algorithm (see IStarck et al.l 1 19981 ; IStarck fc Murtaghl 
2006), which simultaneously computes the wavelet trans- 
form on dyadic scales corresponding to 2 J pixels. This al- 
gorithm decomposes the convergence map of size N x N 
into J = jmax + 1 sub-arrays of size iV x N as follows: 



n(x,y) = cj(x,y) + ^ Wj(x,y) 



(9) 



where j ma x represents the number of wavelet bands (or, 
equivalently, aperture mass maps) considered, cj represents 
a smooth (or continuum) version of the original image k, 
and Wj represents the input map filtered at scale 2 J (i.e. the 
aperture mass map at 6 — 2 J pixels). 

Using the wavelet formalism to derive the aperture mass 
statistic presents different advantages: 

• Many families of wavelet functions have been studied 
in the statistical literature, and all these wavelet functions 
could be applied to the aperture mass statistic. We will 
demonstrate in § U that the starlet wavelet, in particular, 
simultaneously meets all three of the requirements for aper- 
ture mass filters described in § [21 and is therefore ideal for 
weak lensing studies. 

• For some specific wavelet functions, discrete and very 
fast algorithms exist, allowing us to compute a set of wavelet 
scales through the use of a filter bank with a very limited 
number of operations. See lStarck et al.l (|2010f ) for a full re- 
view of the different wavelet transform algorithms. 

• Wavelets have shown to be very powerful f or non Gaus- 



sianity studies for bo t h we ak lensing data ( Starck et al 
200i 



200 



Jin et al 



McE wen et al.1 



Pires et al.ll2009l. |2012|) and CMB data dVielva et al 



2005; Starck et al.1 12004 IVielva et alj 12006: 



20081 ; ICurto et al.ll201ll 1 



4 FILTER FUNCTIONS 

We have shown above that a 2D wavelet transform of a map 
of the convergence n is formally identical to computing aper- 
ture mass maps of the convergence - or, equivalently, the 
shear - at the chosen scales. We now consider the functional 
form of the compensated filter functions commonly associ- 
ated with aperture mass studies, and compare with a wavelet 
transform that has been fre quently used for weak lensing 
studies: th e starlet wavelet (|Starck et al.1 120061 ; IPires et al] 
200^ 120121 ). 
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4.1 Aperture Mass Filter Functions 

I Schneider et al.l (| 199ST > have proposed a family of polynomial 
filter functions defined by: 



QW = 



{l + 2) 2 

n e 2 v~ 
(i+£)(2+e) $ 2 



1 + 2 e 2 
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2 



if 



2\ I 



(10) 



where H(9 — i?) is the Heaviside step function, which takes 

the value of 1 for i? ^ 6 and zero for $ > 8. 

More recently , several authors (e.g. Ivan Waerbekdll99g| ; 
Ijarvis et al.l I2004T ) have advocated a filter function of the 
form: 



A 
2 
A_bf_ 

e 2 ! 2 



exp 



exp 

W 2 
"ff 2 



" e 2 



(ii) 



where various choices for the constant b and the over- 
all normalisation A have been used in the literature (see, 
e.g.. ICrittenden et alj|2002l ; IZhang et all 120031 ; Ijarvis et ail 
|2004|). Of specific i nteres t to this paper are the forms used 
by Ivan Waerbekel (|l99S ), who takes b = 4, and that of 
Jarvi s et all (|2004h and Izhang et~ail (|2Q03h . who take b = 



1/20 



Zhang et all (2003) have shown that the form in equa- 



tion J^2j§_oj>limal for we ak lensing skewness measurements, 
and Ivan Waerbekel (| 1998h notes that this f o rm is preferable 
to the polynomial filter of 1 Schne ider et alJ (|1998T ) for skew- 
ness measurements because the power spectrum of UijD) is 
non-oscillatory, with a well-localised peak in the frequency 
domain. 

Note that the Q(i9) filter function described in equation 
(fTT|) shows a peak at ■& = V20, and tends to zero as # — > oo. 
In practice, any algorithm used to generate an aperture mass 
map will need to truncate this filter function at some finite 
radius. This will involve a trade-off between accuracy and al- 
gorithm speed, and truncation may affect the effective width 
of the filter. This will be discussed more fully in sections l4.3l 
4~4l and [5] below. 



4.2 Starlet Wavelet Function 

One of the most popular wavelet transform algorithms in 
astronom y is the Isotropic Undec imated Wavelet Tranform 
(IUWT) (|Starck et aT]|2010l . l2oTlT l. also called a trous algo- 
rithm or starlet transform. In this algorithm, the wavelet 
ip(x, y) is separable and can be defined by: 



*(f.f)-« 



¥>(z,2/) 



1 fx y\ 
4^2-1) 



(12) 



where f(x, y) — tp(x)ip(y) and <p(x) is a scaling function 
from which the wavelet is generated. In the case of the starlet 



This form is analogous to the Mexican Hat wavelet function 



wavelet, tp(x) is a B3-spline: 

<f(x) = ^(|a;-2| 3 -4| a ;-l| 3 +6|x| 3 -4|x + l| 3 + |x + 2| 3 ), 

(13) 

which is a compact function that is identically zero for \x\ > 
2. 

This wavelet function has a compact support in real 
space, is well localized in Fourier domain, and the wavelet 
decomposition of an image can be obtain ed with a very fast 
algorithm (see IStarck fc Murtaghl (|20061 ) for a full descrip- 
tion). Computation time will be discussed in detail in Sec- 
tion [5] 

Evaluation of equation (|12|l yields the following aper- 
ture mass filters: 
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4.3 Properties in Real Space 

In Figure [T] we compare the aperture mass filter functions 
defined above with the corresponding starlet wavelet func- 
tions in real space. The label S98 corresponds to the filters 
in to equation (|10p taking I = 3, and JBJ04 corresponds to 
equation (|ll|l with b — 1/2. All filters are renormalised to 
have a peak value of 1, for compa rison purposes. We h ave 
not included the filters defined by Ivan Waerbekel (|l998T) 111 
th is comparison, as t hey are functionally identical to those 
arvis et al.l l|2004h with a different scaling. In practice, 
they are very similar in scale and shape to the starlet wavelet 
filters. 

While the S98 t/(i?) filter is compensated within an 
aperture of radius 6, thus yielding an obvious truncation 
point for this filter function when making maps, neither the 
starlet or JBJ04 filters show such behaviour. The starlet fil- 
ter U ($) filter is, however, compensated within an aperture 
of radius 29, and application of the starlet transform algo- 
rithm considers the full filter function out to this radius. In 
contrast, the JBJ04 Q(i9) function tends to zero as 1? — > 00, 
and the filter function is only strictly compensated at 

this radius. 

In practice, a truncation at large e.g. i? CU i ~ 58 would 
yield minimal error on the solution. However, the larger the 
truncation radius, the longer the computation time when 
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Figure 1. Comparison of the U{9/6) (left) and Q{$/6) (right) 
filter functions commonly used in weak lensing studies and corre- 
sponding starlet wavelet functions. Plots are normalised so that 
each curve attains a peak amplitude of 1. The filters considered 
are those of Schneider et al. (1998) (taking I = 3, labelled S98), 
and Jarvis et al. (2004) (equation (JTTJ with b = 1/2, labelled 
JBJ04). A dotted vertical line in each plot shows the location of 
the aperture radius, indicating where the filters would be trun- 
cated in application of the aperture mass technique. 



computing an aperture mass map. Moreover, contamination 
by edge effects and missing data becomes more problematic 
as the aperture radius is increased. Two-point and higher 
order statistics, however, remain unaffected by this trunca- 
tion, as the effect of the aperture mass filter can be com- 
puted analytically in both real and Fourier space for these 
statistics. 

Moreover, statistics such as peak counting are unlikely 



to be affected by such a truncation. For a non-compensated 
filter function, the amplitude of structures seen in the result- 
ing Map realisation may be offset by an additive constant 
proportional to /*«' §U(§)dd (see ISchneid"eilll996t ). Peak 
counting is generally performed on the signal to noise of the 
M a p, and as the noise map is generated in the same manner 
as the signal, it stands to reason that this will not be af- 
fected significantly. A similar argument applies to statistical 
moments of the M ap distribution such as the variance, skew- 
ness and kurtosis computations performed on the aperture 
mass map directly, which subtract the mean from the map 
inherently. 



4.4 Properties in Fourier Space 

Perhaps more important than their real-space behaviour is 
the behaviour of these filter functions in Fourier space. The 
key feature of aperture mass methods is to isolate structures 
on a particular scale determined by the aperture radius. If 
this radius is chosen to be much larger than the scale at 
which the noise is dominant, one may suppress the noise and 
therefore resolve real structures above the noise. In this way, 
the aperture mass filters (or, equivalently, the wavelet trans- 
form) act as band-pass filters, restricting a given aperture 
mass map or wavelet scale to a limited range of frequencies 
in the Fourier domain. 

To assess the response of our aperture mass filters and 
starlet wavelet transform in Fourier space, we consider the 
power spectrum of the U(§) filter function. Ideally, we would 
like the chosen filter to be localised in the Fourier domain, 
with a single peak corresponding to the characteristic fre- 
quency for the scale being considered, and with no oscil- 
lations at high frequencies, as these would result in high- 
frequency (noisy) information being included in a lower- 
frequency filtered map. 

In order to accurately characterise the behaviour of the 
filters, including any truncation that is applied, we generate 
an artificial shear data generated from a null convergence 
map with a single central delta function, and with no noise 
in the maps. The aperture mass algorithm is then applied 
to the resulting shear maps for 9 — 4, and the power spec- 
trum of the resulting map computed. For the JBJ04 filter, 
we consider two truncation radii, at ■dcut = 9 and $ cut = 58, 
for illustrative purposes. The wavelet transform is also com- 
puted for the map, and the power spectrum for the second 
wavelet scale computed for comparison. The power spectra 
are then renormalised to peak at 1. Figure [2] shows the re- 
sulting power spectra. In the figure, each curve has been 
offset from the one below by 0.5 along the vertical axis, for 
clarity. 

Broadly, all the filter functions show the expected be- 
haviour. All but one are peaked at approximately the same 
frequency, indicative of the fact that they are probing struc- 
tures at a similar scale, and they are all fairly well localised 
in the frequency domain. The frequency offset seen in the 
peak of the JBJ04 filter truncated at ■& = 58 is indicative 
of the change in the effective width of the filter function as 
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Figure 2. Power spectra associated with the aperture mass filters 
S98, JBJ04 and the starlet wavelet filter. We have plotted the S98 
for both £ = 3 as before (solid curve) and 1 = 1 (dotted curve) to 
show the dependence of the power spectrum on the polynomial 
order I. For the JBJ04 filter, a truncation has been applied at 
$<m< = 58 (solid curve) and i? ctl t = (dot-dashed curve). The 
power spectra for the different filter functions have been offset at 
intervals of 0.5 along the vertical axis, for clarity. 



the truncation radius is changed. Oscillations are seen with 
the JBJ04 filter truncated at $ cut = 6 and the S98 filter at 
high frequency, but not in the case of the starlet filter or the 
JBJ04 filter truncated at i? cut = 59. 

The oscillations in the strongly truncated JBJ04 case 
are a result of the fact that the filter function is significantly 
non-zero and non-compensated within the truncation radius 
of &cut = 0. This results in oscillations at high frequency 
and, as a consequence, a contamination of the resulting M ap 
map with high frequency modes not pertaining to the scale 
of interest. Clearly, truncations at a large radius, such as 
$cut = 58 yield negligible oscillations. 

Similar oscillations are seen in the S98 filter. This is a 
direct consequence of the fact that the filter used is a locally- 
defined function: it is defined to follow the given form only 
for § < 6 and artificially set to zero outside this radius by 
the Heaviside step function. This step function introduces 
a discontinuity in the filter function, which manifests as os- 
cillations in the Fourier space response of the filter. These 
oscillations are reduced in amplitude as I is increased. 

The starlet filter is very well-behaved at high frequen- 
cies, showing no oscillations at any frequency, which results 
from the fact that the starlet transform algorithm applies no 
truncation to the starlet filter functions. Therefore, the star- 
let transform at each scale represents a strictly limited range 
of frequencies, with no contamination from higher- frequency 
modes. 

The starlet filter offers an additional advantage over 
aperture mass filtering: the aperture mass filter is a band- 
pass filter. Because the M ap reconstruction is only computed 
for a discrete set of aperture scales, information on interme- 
diate scales may be lost, particularly at both the small and 
large frequency extremes. In contrast, the wavelet transform 
retains information on all scales. The first wavelet scale is 
effectively a high-pass filter, retaining all the high-frequency 
information in the image, while the remaining wavelet scales 
are bandpass filtered as with the M ap -filtered images. Fi- 
nally, the wavelet transform retains cj, which encodes the 
large-scale information in the image, and therefore consists 
of a smoothed version of the image (or an image with a low- 
pass filter applied). Figure [3] demonstrates this point for a 
wavelet transform with J = 5, showing the Fourier-space 
response of the wavelet filter as a function of scale j. 

4.5 Practical considerations in application 

The starlet transform would seem to offer a clear advantage 
over traditional M ap filtering techniques when it comes to 
map making, as it requires no truncation of the filter func- 
tions and retains information on all scales within the original 
image. However, in application there are several difficulties 
that arise, and which must be addressed. 

The starlet transform is computed on the convergence 
map, rather than on the shear catalogue itself. The shear 
and convergence are related by a straightforward 2D convo- 
lution which, in principle, may be inverted directly in Fourier 
space. However, in practice, this gives rise to significant er- 
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Figure 3. Fourier-space response of the wavelet filter as a func- 
tion of filter scale. Note that at at intermediate scales, the wavelet 
filter responds as a limited band-pass filter, while at large and 
small scales, the wavelet transform yields a low- and high-pass 
filter, respectively. 



rors in the convergence estimation, particularly at the edges 
of the field, arising from the fact that we are considering 
a finite field sampled over a discrete range of frequencies, 
while the convolution acts over all space. 

This can give rise to ~ 10% leakage of power into lensing 
B-modes in the convergence map, which biases the resulting 
convergence map. Various techniques have been developed 
to mitigate such effects; most recently, iDeriaz et ail (|2012h 
have presented a method using a Wavelet-Helmholtz decom- 
position to separate E- and B-modes in the lensing signal, 
and to recover the convergence to excellent precision, of- 
fering a factor of ~ 30 x reduction in B-mode contamina- 
tion compared to a direct FFT inversion and ~ 5x reduc- 
tion in B-mode contamination compared to the method of 
ISeitz fc Schneider! ljl996l . l200ll ). and yielding an RMS error 
at the percent leve l in reconstructions. Note also that the 
IDeriaz et al.l l|2012r ) method may be applied directly to dis- 
crete data; i.e. the technique can be applied directly to a 
shear catalogue, rather than pixelated shear maps. 

A further difficulty that arises when attempting to re- 
construct the convergence, and compute statistics of the 
resulting map, is due to missing data in the image. This 
usually arises due to masking out of bright stars and other 
contaminating features in the image (such as bad pixels). 
This lost information is not recovered in the reconstruction, 
and this biases the statistic s computed on the resulting mass 
map. In iPires et all |2009r i, the authors describe a method 



by which one can recover information in regions covered by 
the mask. This method is based on an assumption of spar- 
sity, and the authors find that they are able to recover the 
power spectrum to a relative error of 0.5 — 1% using Sub- 
aru and CFHT masks, respectively on a 4 deg 2 image, and 
can recover the equilateral bispectrum to a relative error of 
1 — 3% for the same masks. 

Using the techniques described above, it is therefore 
possible to compute maps of the convergence to sufficient 
accuracy to be able to compute associated statistics on the 
resulting maps, and thereby place constraints on cosmolog- 
ical parameters in a similar manner to methods applied di- 
rectly to the shear catalogue. 



5 ALGORITHM SPEEDS 
5.1 Map making 

Generating an aperture mass map generally involves brute 
force application of equation |2[ to a given dataset. Such 
an algorithm has complexity oc 0{N 2 fl 2 ut ), where N x TV is 
the dimension of the image, and # cu t is the chosen trunca- 
tion radius. This scaling means that for large apertures or, 
equivalently, for high-resolution images, the algorithm may 
prove to be very time-consuming. Estimation of the noise in 
aperture mass studies is generally carried out by randomis- 
ing the data and repeating the measurement many times. 
Therefore it is of great importance to be able to compute 
the aperture mass for a given image rapidly. 

The starlet wavelet transform algorithm is of complexity 
oc 0(N 2 J), where J is the number of scales considered, and 
is limited by N ^ 2 J . This means that the processing time 
for the wavelet transform algorithm is sensitive only to the 
number of scales considered, rather than the size of the filter 
functions involved, and depends linearly on this number. 

In Figure 21 we compare the processing time for the 
aperture mass algorithm and the starlet transform algo- 
rithrrjf], both programmed in CH — h, to analyse an image of 
1024 x 1024 pixels on a 2 x 2.66GHz Intel Xeon Dual-Core 
processor. We consider aperture scales 9 = [2, 4, 8, 16, 32, 64], 
which correspond to J = j max + 1 = [2, 3, 4, 5, 6, 7] wavelet 
scales in the wavelet transform. In the aperture mass al- 
gorithm, the filters are truncated at a radius of _> cu t = 9. 
For the application of filters such as the JBJ04 filter, which 
necessitates truncation at a much larger radius, we expect 
the computation time to be roughly an order of magnitude 
longer. Even at the smallest aperture radius, the wavelet 
transform is ~ 5x faster than the aperture mass algorithm. 
At 9 — 64 pixels, the wavelet transform is ~ 1200 x faster 
than the aperture mass algorithm. Note that the wavelet 
transform for a given value of J simultaneously computes 
the wavelet transform at all scales 2 3 , < j ^ J — 1, in addi- 
tion to the smoothed continuum map cj, whilst the aperture 
mass algorithm computes the transform at a single scale 9. 

3 available from http://irfu.cea.fr/Ast/mrlens_software.php 
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Figure 4. Comparison of the processing time of the aperture 
mass algorithm and the starlet transform algorithm to analyse 
an image of 1024 X 1024 pixels on a 2 X 2.66GHz Intel Xeon 
Dual-Core processor. The top axis represents the total number J 
of 1024 X 1024 pixel arrays computed by the wavelet transform 
algorithm, including the smooth continuum map. 



We note further that the computational time for the wavelet 
transform for J = 7 wavelet scales is still a factor of ~ 2x 
less than the computational time for the aperture mass al- 
gorithm at 6 = 2 pixels. 

The starlet transform algorithm therefore offers a clear 
and significant time advantage over the aperture mass al- 
gorithm for all scales of interest. Note that wavelets do not 
restrict our data analysis only to dyadic scales. Indeed some 
other wavelet tran sform algorithms, such as that based on 
fractional splines (|Unser &: Blul l200(jh . allow us to use any 
intermediate wavelet scale. 



5.2 Computing higher-order M ap statistics 

When computing higher-order statistics such as the vari- 
ance, skewness or kurtosis of the aperture mass statistic, one 
can choose to work directly on the shear field, and compute 
the statistics via the n-point correlation function. Computa- 
tion of such statistics using a naive approach usually results 
in an algorithm with complexity 0(N" ;), where N ga i is the 
number of galaxies being considered and n is the order of 
the correlation function. On large fields, such a computation 
can be prohibitive. 

In recent years, tree codes have been employed to speed 
up computation of n-point correlation functions. Typical 
tree codes to com pute n-point correl ation functions are 
0(N ga i log(N gal )) \ Zhang fc Pen! 120051 ) on a shear cata- 
logue. For a single Euclid exposure of 0.5 deg 2 , we can ex- 
pect Ng a i ~ 54,000 (30 galaxies/arcmin 2 ) - 180,000 (100 
galaxies/a rcmin 2 ). Tree codes exist that act on pixelated 
data (e.g. Eri ksen et al.l 2004) which run at 0(Np ix n oin ) 
where Np ix is the total number of pixels and nun is the num- 
ber of bins in the correlation function. For a Euclid exposure, 



assuming pixels of 1 arcminute, we have N pix = 1800, and 
ribin will be dependent on the required resolution of the cor- 
relation function. Typically, statistics wi ll be computed on 
much larger fields, however (for example. IPires et al.l |2009l ) 
compute statistics on a field of 3.95 x 3.95 deg 2 , which is a 
scale- up of ~ 30 in area). 

The wavelet method acts on pixellated data, and is 
0(N 2 J) in computation time, so our algorithm will be com- 
parable for computation of 2-point statistics, if the 2-point 
correlation function is computed on pixellated data, but if 
a shear catalogue is used, we will have a faster algorithm 
by at least an order of magnitude. For higher-order statis- 
tics, this advantage is even more pronounced. Furthermore, 
while optimised software is freely and publicly available to 
compute the wavelet transform, such optimised software is 
not available f or n-point cor r elatio n functions. 

Note that Ijarvis et all (2004) show an 0{N ga i) com- 
plexity for the computation of n-point statistics with their 
algorithm, with no log(iV ga ;) factor. However, for n ^ 3, the 
prefactor becomes quite large, and for n = 3 their algorithm 
was unable to reach the theoretical 0(N ga i log(N ga i)) scal- 
ing at N ga i = 10 6 . This means that their algorithm may be 
competitive for computation of two-point statistics, but for 
higher order statistics, methods acting on pixellated data 
- such as the wavelet method described above - would be 
advantageous. 



5.3 Further advantages of wavelets 

In lStarck et al.l |2006l ). the authors describe a method for re- 
constructing the convergence from shear measurements us- 
ing wavelets. Using their method, it is possible to carry out 
an explicit denoising of the convergence field using thresh- 
olding based on a False Discovery Rate method. This allows 
them to derive robust detection levels in wavelet space, and 
to produce high-fidelity denoised mass maps. In lPires et al.l 
(|2009l . 120121 ), it is shown that more non-gaussian informa- 
tion can be extracted from these wavelet-denoised maps as 
compared to the shear, convergence or M ap -filtered maps, 
and therefore that tighter constraints may be placed on the 
cosmological model using the wavelet method. 

In addition, wavelet-based methods offer more flexibil- 
ity than aperture mass filters. Whilst we have presented here 
only the starlet wavelet function, many other wavelet dictio- 
naries may be used. The starlet filter seems ideal for lensing 
studies, due to its similarity to the J BJ04 filter presente d 
here, which was deemed to be optimal in lZhang et al.l (2003). 
However, different dictionaries may be optimal in different 
applications; for example, if one were attempting to study 
filamentary structure, ridgelets or curvelets might be a more 
appropriate basis. The vastness of the wavelet libraries and 
the public availability of fast algorithms to compute these 
transforms are major strengths of wavelet-based approaches. 
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6 SUMMARY AND DISCUSSION 

In this paper, we have compared the aperture mass statistic, 
a commonly-used tool in weak lensing studies, with the 2D 
wavelet transform using a starlet wavelet. We have demon- 
strated that the aperture mass technique is formally identi- 
cal to the wavelet transform evaluated at a specific scale. 

We have considered two common forms for the aperture 
mass filter function, and compared it to the equivalent star- 
let functions. One of these aperture mass filters - labelled 
JBJ04 - is considere d to be the optimal aperture mass filter 
for skewness studies l|Zhang et al.| [2003). and is thus consid- 
ered to be the state-of-the-art in this field. We outline three 
criteria for the aperture mass filters: 

(i) They should be local in real space. 

(ii) They should be compensated within the aperture ra- 
dius being considered. 

(iii) They should be local and non-oscillatory in Fourier 
space. 

We demonstrate that none of the commonly used aperture 
mass filter functions meet all three criteria simultaneously, 
whilst the starlet - by definition - meets all of the criteria. 

Lastly, we show that the starlet transform algorithm, 
which forms part of a publicly available software package 
containing numerous tools for weak lensing, dramatically 
outperforms a brute force application of the aperture mass 
technique, with speed-up factors of 5— ~ 1200 seen for the 
processing of a single aperture mass reconstruction, depend- 
ing on the scale considered. 

Given that the aperture mass technique is identical to 
the wavelet transform, that the starlet wavelet function is 
more appropriately suited to aperture mass studies, and that 
the starlet transform algorithm offers such a great advantage 
in terms of processing time over the standard aperture mass 
algorithm, we argue that the wavelet transform should be 
the preferred method in future weak lensing studies involv- 
ing use of the aperture mass method. 

A further advantage of this wavelet method is that 
fast and effective software is availa ble to carry out exp licit 
denoising in the wave let domain l|Starck et al.l I2006T ). In 
iPires et al. I (|2009l , 120121 ) . it has been shown that further cos- 
mological information can be extracted from noisy conver- 
gence maps if such denoising is applied to the convergence 
maps, as the starlet wavelet transform is particularly sen- 
sitive to non-gaussian structures. In addition, a multitude 
of wavelet libraries are available, and may be chosen to be 
optimal for a given problem. Thus, wavelet methods have 
an inherent flexibility in application. 
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